Sensitivity and Generalization of a Neural Network for Estimating Left Atrial Fibrotic Volume Fractions from the 12-lead ECG

Features extracted from P waves of the 12-lead electrocardiogram (ECG) have proven valuable for non-invasively estimating the left atrial fibrotic volume fraction associated with the arrhythmogenesis of atrial fibrillation. However, feature extraction in the clinical context is prone to errors and oftentimes yields unreliable results in the presence of noise. This leads to inaccurate input values provided to machine learning algorithms tailored at estimating the amount of atrial fibrosis with clinical ECGs.Another important aspect for clinical translation is the network’s generalization ability regarding newECGs.To quantify a network’s sensitivity to inaccurately extracted P wave features, we added Gaussian noise to the features extracted from 540,000 simulated ECGs consisting of P wave duration, dispersion, terminal force in lead V1, peak-to-peak amplitudes, and additionallythoracic and atrial volumes. For assessing generalization, we evaluated the network performance for train-validation-test splits divided such that ECGs simulated with the same atria or torso geometry only belongedto either the trainingand validationor the test set. The root mean squared error (RMSE) of the network increased the most in case of noisy torso volumes and P wave durations. Large generalization errors witha RMSEdifference between training and test set of more than 2% fibrotic volume fraction only occurred ifveryhigh or low atria and torso volumes were left out during training.Our results suggest that P wave duration and thoracic volume are features that have to be measured accurately if employed for estimating atrial fibrosis with a neural network. Furthermore, our method is capable of generalizing wellto ECGs simulated with anatomical models excluded during training and thus meets an important requirement for clinical translation

Electrocardiographic Imaging Using a Spatio-Temporal Basis of Body Surface Potentials - Application to Atrial Ectopic Activity

Electrocardiographic imaging (ECGI) strongly relies on a priori assumptions and additional information to overcome ill-posedness. The major challenge of obtaining good reconstructions consists in finding ways to add information that effectively restricts the solution space without violating properties of the sought solution. In this work, we attempt to address this problem by constructing a spatio-temporal basis of body surface potentials (BSP) from simulations of many focal excitations. Measured BSPs are projected onto this basis and reconstructions are expressed as linear combinations of corresponding transmembrane voltage (TMV) basis vectors. The novel method was applied to simulations of 100 atrial ectopic foci with three different conduction velocities. Three signal-to-noise ratios (SNR) and bases of six different temporal lengths were considered. Reconstruction quality was evaluated using the spatial correlation coefficient of TMVs as well as estimated local activation times (LAT). The focus localization error was assessed by computing the geodesic distance between true and reconstructed foci. Compared with an optimally parameterized Tikhonov-Greensite method, the BSP basis reconstruction increased the mean TMV correlation by up to 22, 24, and 32% for an SNR of 40, 20, and 0 dB, respectively. Mean LAT correlation could be improved by up to 5, 7, and 19% for the three SNRs. For 0 dB, the average localization error could be halved from 15.8 to 7.9 mm. For the largest basis length, the localization error was always below 34 mm. In conclusion, the new method improved reconstructions of atrial ectopic activity especially for low SNRs. Localization of ectopic foci turned out to be more robust and more accurate. Preliminary experiments indicate that the basis generalizes to some extent from the training data and may even be applied for reconstruction of non-ectopic activity

A bi-atrial statistical shape model for large-scale in silico studies of human atria: model development and application to ECG simulations

Large-scale electrophysiological simulations to obtain electrocardiograms (ECG) carry the potential to produce extensive datasets for training of machine learning classifiers to, e.g., discriminate between different cardiac pathologies. The adoption of simulations for these purposes is limited due to a lack of ready-to-use models covering atrial anatomical variability. We built a bi-atrial statistical shape model (SSM) of the endocardial wall based on 47 segmented human CT and MRI datasets using Gaussian process morphable models. Generalization, specificity, and compactness metrics were evaluated. The SSM was applied to simulate atrial ECGs in 100 random volumetric instances. The first eigenmode of our SSM reflects a change of the total volume of both atria, the second the asymmetry between left vs. right atrial volume, the third a change in the prominence of the atrial appendages. The SSM is capable of generalizing well to unseen geometries and 95% of the total shape variance is covered by its first 23 eigenvectors. The P waves in the 12-lead ECG of 100 random instances showed a duration of 104ms in accordance with large cohort studies. The novel bi-atrial SSM itself as well as 100 exemplary instances with rule-based augmentation of atrial wall thickness, fiber orientation, inter-atrial bridges and tags for anatomical structures have been made publicly available. The novel, openly available bi-atrial SSM can in future be employed to generate large sets of realistic atrial geometries as a basis for in silico big data approaches

A Reproducible Protocol to Assess Arrhythmia Vulnerability in Silico: Pacing at the End of the Effective Refractory Period

In both clinical and computational studies, different pacing protocols are used to induce arrhythmia and non-inducibility is often considered as the endpoint of treatment. The need for a standardized methodology is urgent since the choice of the protocol used to induce arrhythmia could lead to contrasting results, e.g., in assessing atrial fibrillation (AF) vulnerabilty. Therefore, we propose a novel method—pacing at the end of the effective refractory period (PEERP)—and compare it to state-of-the-art protocols, such as phase singularity distribution (PSD) and rapid pacing (RP) in a computational study. All methods were tested by pacing from evenly distributed endocardial points at 1 cm inter-point distance in two bi-atrial geometries. Seven different atrial models were implemented: five cases without specific AF-induced remodeling but with decreasing global conduction velocity and two persistent AF cases with an increasing amount of fibrosis resembling different substrate remodeling stages. Compared with PSD and RP, PEERP induced a larger variety of arrhythmia complexity requiring, on average, only 2.7 extra-stimuli and 3 s of simulation time to initiate reentry. Moreover, PEERP and PSD were the protocols which unveiled a larger number of areas vulnerable to sustain stable long living reentries compared to RP. Finally, PEERP can foster standardization and reproducibility, since, in contrast to the other protocols, it is a parameter-free method. Furthermore, we discuss its clinical applicability. We conclude that the choice of the inducing protocol has an influence on both initiation and maintenance of AF and we propose and provide PEERP as a reproducible method to assess arrhythmia vulnerability

Separating ventricular activity in thoracic EIT using 4D image-based FEM simulations

One challenge in central hemodynamic monitoring based on electrical impedance tomography (EIT) is to robustly detect ventricular signal components and the corresponding EIT image region without external monitoring information. Current stimulation and voltage measurement of EIT were simulated with finite element porcine torso models in presence of a multitude ofthoracic blood volume shifts. The simulated measurement data was examined for linear dependence on changes in stroke volume. Based onthe results the EIT measurement information regardingstroke volume changesis sparse

Refining the Eikonal Model to Reproduce the Influence of Atrial Tissue Geometry on Conduction Velocity

Atrial fibrillation is responsible for a significant and steadily rising burden. Simultaneously, the treatment options for atrial fibrillation are far from optimal. Personalized simulations of cardiac electrophysiology could assist clinicians in the risk stratification and therapy planning for atrial fibrillation. However, the use of personalized simulations in clinics is currently not possible due to either too high computational costs or non-sufficient accuracy. Eikonal simulations come with low computational costs but cannot replicate the influence of cardiac tissue geometry on the conduction velocity of the wave propagation. Consequently, they currently lack the required accuracy to be applied in clinics. Biophysically detailed simulations on the other hand are accurate but associated with too high computational costs. To tackle this issue, a regression model is created based on biophysically detailed bidomain simulation data. This regression formula calculates the conduction velocity dependent on the thickness and curvature of the heart wall. Afterwards the formula was implemented into the eikonal model with the goal to increase the accuracy of the eikonal model without losing its advantage of computational efficiency. The results of the modified eikonal simulations demonstrate that (i) the local activation times become significantly closer to those of the biophysically detailed bidomain simulations, (ii) the advantage of the eikonal model of a low sensitivity to the resolution of the mesh was reduced further, and (iii) the unrealistic occurrence of endo-epicardial dissociation in simulations was remedied. The results suggest that the accuracy of the eikonal model was significantly increased. At the same time, the additional computational costs caused by the implementation of the regression formula are neglectable. In conclusion, a successful step towards a more accurate and fast computational model of cardiac electrophysiology was achieved